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Mean-field theory predicts that bilayer quantum Hall systems at odd integer total filling factors 
can have stripe ground states in which the top Landau level is occupied alternately by electrons 
in one of the two layers. We report on an analysis of the properties of these states based on a 
coupled Luttinger liquid description that is able to account for quantum fluctuations of charge- 
density and position along each stripe edge. The soft modes associated with the broken symmetries 
of the stripe state lead to an unusual coupled Luttinger liquid system with strongly enhanced low- 
temperature heat capacity and strongly suppressed low-energy tunneling density of states. We assess 
the importance of the intralayer and interlayer back-scattering terms in the microscopic Hamiltonian, 
which are absent in the Luttinger liquid description, by employing a perturbative renormalization 
group approach which re-scales time and length along but not transverse to the stripes. With 
interlayer back-scattering interactions present the Luttinger liquid states are unstable either to an 
incompressible striped state that has spontaneous interlayer phase coherence and a sizable charge 
gap even at relatively large layer separations, or to Wigner crystal states. Our quantitative estimates 
of the gaps produced by back-scattering interactions are summarized in Fig. 11 by a schematic phase 
diagram intended to represent predicted experimental findings in very high mobility bilayer systems 
at dilution refrigerator temperatures as a function of layer separation and bilayer density balance. 
We predict that the bilayer will form incompressible isotropic interlayer phase coherent states for 
small layer separations, say d < 1.5£. At larger interlayer spacings, however, the bilayer will tend 
to form one of several different anisotropic states depending on the layer charge balance, which we 
parameterize by the fractional filling factor v contributed by one of the two layers. For large charge 
imbalances (y far from 1/2), we predict states in which anisotropic Wigner crystals form in each 
of the layers. For v closer to 1/2, we predict stripe states that have spontaneous inter-layer phase 
coherence and a gap for charged excitations. These states should exhibit the quantum Hall effect for 
current flowing within the layers and also the giant interlayer tunneling conductance anomalies at 
low bias voltages that have been observed in bilayers when the N = Landau level is partially filled. 
When the gaps produced by backscattering interactions are sufficiently small, the phenomenology 
observed at typical dilution fridge temperatures will be that of a smectic metal, anisotropic transport 
without a quantum Hall effect. For stripe states in the N — 2 Landau level, this behavior is expected 
over a range of bilayer charge imbalances on both sides oiv= 1/2. 



I. INTRODUCTION 

The recent discovery of strongly anisotropic transport 
in single layer quantum Hall systems near half-odd in- 
teger filling factors 1-3 has attracted much experimental 4 
and theoretical 5 interest. Transport anisotropies have 
been observed in single two-dimensional (2D) electron gas 
layers at half filling of Landau levels with index N > 2, 
i.e. at filling factors v = 9/2,11/2,.... This effect is 
commonly ascribed to the formation of striped charge- 
density-wave phases, predicted on the basis of Hartree- 
Fock calculations by Koulakov et al. e and by Moessner 
and Chalker 7 with additional theoretical support from 
subsequent exact-diagonalization 8 and DMRG 9 numer- 
ical studies. The stripe state is a consequence of the 
form factors that arise in describing interactions be- 
tween electrons in higher kinetic energy Landau level 
orbitals and allow density waves to form in cyclotron- 
orbit-center coordinates that have a very small electron- 



density-wave amplitude and therefore little electrostatic 
energy penalty. 

The physics of quantum Hall systems is enriched by 
the additional degrees of freedom that appear in bilayer 
systems 10 in which two 2D electron layers have a sepa- 
ration d small enough that their interactions have con- 
sequences. For total filling factor vt = 1 and other 
odd integer total filling factors, interlayer interactions 
can lead to a state with spontaneous phase coherence 11 
between the layers and a charge gap that is revealed 
experimentally 12 by the quantum Hall effect. Further 
spectacular experimental manifestations of spontaneous 
phase coherence were revealed very recently in 2D to 2D 
tunneling and Hall drag experiments by Eisenstein and 
collaborators. 13 . In tunneling studies spontaneous coher- 
ence is signaled by a sharp zero bias peak in the differ- 
ential conductance between the layers. As the ratio of d 
to the magnetic length i is reduced experimentally, the 
conductance peak appears to develop continuously start- 



1 



ing at a critical value of djl that is consistent with ear- 
lier experimental anomalies 12 attributed to spontaneous 
coherence and with mean-field-theory estimates of the 
critical layer separation 12 at which coherence is expected 
to develop. These experiments are still not understood 
quantitatively and raise a number of interesting issues 
in non-equilibrium collective transport theory that have 
stimulated a growing body of theoretical 14 work. 

Since balanced bilayer systems at large odd integer to- 
tal filling factor (i>t > 9) are composed of 2D layers that, 
if isolated, would show stripe-state behavior, it is natu- 
ral to consider the possible interplay and competition 
between the formation of striped phases in each 2D layer 
and the development of spontaneous interlayer phase- 
coherence. These issues have been investigated in several 
recent theoretical papers 15-17 and it has been argued 17 
that they may be relevant for understanding a recent ob- 
servation of resistance anisotropy at integer filling factor 
by Pan et a/. 18 . In the present paper we extend earlier 
work by two of the present authors 19 on smectic states in 
single layers to the case of bilayer systems. The approach 
we take is one that is intended to be valid when quan- 
tum fluctuation corrections to the stripe states predicted 
by Hartree-Fock theory are weak on microscopic length 
scales, although as we discuss at length they inevitably 
alter the ultimate physics at very low energies and tem- 
peratures and the behavior of correlation functions at 
long distances. Since stripe states occur as extrema of 
the Hartree-Fock energy functional for any orbital Lan- 
dau level index, not only for N > 2 where the states 
are seen experimentally, and are in fact always unsta- 
ble to the formation of Wigner crystal states in mean- 
field-theory, it is evident that we must appeal in part 
to experiment to judge when our starting assumption is 
valid. 5 - 6 ' 20 

In describing stripe states it is convenient to use a 
Landau gauge basis with single-particle states extended 
in the direction along the stripes, which wc choose to 
be the x direction, and labeled by a one-dimensional 
wavevector k that is proportional to the guiding center 
along which the wavefunction's y-coordinate is localized, 
Y = k£ 2 . For balanced bilayers, the stripe states that oc- 
cur in Hartree-Fock theory are occupation number eigen- 
states in this representation, with the valence Landau- 
level Landau gauge states occupied by top and bottom 
layer electrons in alternating stripes. In the Hartree- 
Fock approximation, the low-energy excitations of the 
stripe states consist of coupled particle-hole excitations 
along each edge of top and bottom layer stripes. These 
degrees of freedom are conveniently described using the 
bosonization techniques familiar from the theory of one- 
dimensional electron systems. 21 Our approach is partly in 
the spirit of Fermi liquid theory in that we assume that 
the Hilbert space of low-energy excitations can be placed 
in one-to-one correspondence with those that occur in 
the Hartree-Fock theory. When quantum fluctuations 
are too strong our approach will not be useful; for exam- 
ple, it cannot predict either the fact that lowest Landau 



level isolated layers have composite-fcrmion liquid rather 
than stripe ground states, or the likelihood of bubble 6 
rather than stripe states far away from half-filling. Our 
approach to stripe state physics is similar to that taken 
first by Fradkin and Kivelson 22 . For the case of monolay- 
ers, the microscopic basis of the coupled Luttinger liquid 
model for quantum Hall stripe states was carefully ex- 
amined by Lopatnikova et al. 23 and other properties of 
quantum Hall stripe states have been addressed by Barci 
et al. 24 and Wexler and Dorsey. 25 

Our paper is organized as follows. In Section II we 
review the coupled Luttinger liquid model for quantum 
Hall stripe states and discuss its application to the bilayer 
case. The model rests fundamentally on the assumption 
that the excitation spectrum of bilayer stripe states may 
be placed in one-to-one correspondence with that of the 
Hartree-Fock picture; for bilayers this assumption im- 
plies that the degrees of freedom at each stripe edge are 
those of a one-dimensional electron gas. Our analysis 
of the low-energy long-wavelength physics examines this 
subspace of the microscopic many-particle Hilbert space 
and includes forward scattering terms in the Hamilto- 
nian that create and destroy particle-hole excitations at 
the stripe edges, and back scattering terms that scat- 
ter electrons between chiral one-dimensional electron- 
gas branches. Since the microscopic amplitude of back- 
scattering processes is weak, they can often be neglected 
at experimentally accessible temperatures. When only 
forward scattering terms are included, the Hamiltonian 
can be solved exactly using bosonization and is formally 
equivalent to that of a system of coupled one-dimensional 
electron gases. The quantum smectic broken symmetry 
character of the electronic state is reflected, however, in 
the coupled Luttinger liquid interaction parameters and 
results in enhanced fluctuations. The properties of this 
bilayer smectic state are discussed in Section III. The be- 
havior of the one-particle Green's functions at the smectic 
fixed point, carefully addressed by Lopatnikova et al. 23 
for the single layer case, is discussed for the bilayer. We 
find that, as in the single-layer case, the one-particle 
Green's function does not exhibit the power law behav- 
ior that is generic for weakly coupled Luttinger liquids 
and instead vanishes faster than any power law at large 
distances, implying strongly suppressed tunneling at low 
energies. The enhanced importance of quantum fluctu- 
ations is a consequence of the invariance of the model's 
Luttinger liquid Hamiltonian under a simultaneous trans- 
lation of all stripes. Back-scattering interactions are ad- 
dressed in Section IV, using a perturbative RG approach. 
As in the single-layer case we find that back-scattering 
interactions are always relevant. The gapless Hartree- 
Fock smectic state cannot be the true ground state in ei- 
ther single-layer or bilayer quantum Hall systems. In- 
stead we conclude that except at relatively large layer 
separations, interlayer interactions induce a ground state 
that has spontaneous interlayer phase coherence. This 
state would be signaled experimentally by the simulta- 
neous occurrence of an integer quantum Hall effect and 
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anisotropic finite-temperature transport, something that 
has not been seen in single-layer systems. Where in- 
tralayer interactions are more important, they drive the 
system to a state with an anisotropic Wigner crystal in 
each layer. We argue that both types of interactions lead 
to charge gaps and to integer quantum Hall effects and 
estimate the size of the resulting energy gaps. According 
to our estimates, the gap created by inter-layer back- 
scattering will be large enough to be observable out to 
surprisingly large layer separations. The effect of finite 
tunneling between the layers is also addressed in Section 
IV. Finally in Section V we discuss several interesting 
theoretical issues that arise from this work. We com- 
ment explicitly on inconsistencies between the conclu- 
sions that have been reached by different researchers on 
the question of smectic state stability in the single-layer 
case. We also address the suggestion 24 that the enhanced 
quantum fluctuations that follow from the broken trans- 
lational symmetry of the starting Hartree-Fock state may 
invalidate our perturbative renormalization group analy- 
sis for back-scattering interactions. 

II. THE MODEL 

A. Coupled Luttinger Liquid Model Energy 
Functional 

In Hartree-Fock theory the smectic bilayer state at to- 
tal filling factor vt = 1 is a single Slater determinant 
where the occupation of guiding-center modes in a Lan- 
dau level of index N € {0, 1,2,.. .} alternates between 
the layers with period a, as depicted schematically in 
Fig. 1. Lower index Landau levels are assumed to be 
frozen in filled states and higher index in empty states, 
allowing them to be neglected in the following. Each 
stripe has chiral left-moving and right-moving branches 
of quasiparticle edge states, localized in opposite layers, 
allowing the low-energy degrees of freedom of each elec- 
tron stripe to be mapped to those of a one-dimensional 
electron gas. We consider the general case of a biased 
double layer system where the width of the stripes in 
one layer is av while that in the other layer is a(l — u), 
v £ [0, 1]. As indicated in the figure, for v ^ 0.5 the sys- 
tem has two types of stripe edges, distinguished by the 
direction of their closest neighbor. In the following we 
refer to a pair of chiral stripe edges one above the other 
in each layer as a rung, and the two closest such rungs as 
a rung pair. 

Small fluctuations in the positions and shapes of the 
stripes can be described in terms of particle-hole excita- 
tions near the stripe edges. The residual interactions, 
ignored in Hartree-Fock theory, which act on energy 
states fall into two classes: forward scattering interac- 
tions which conserve the number of electrons on each 
edge of every stripe, and backward scattering processes 
which do not. The quantum smectic model includes 
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FIG. 1. Schematic illustration of the Hartrcc-Fock bilayer 
smectic state. The shaded areas are electron stripes whose 
edges are chiral Luttinger liquids as denoted by the solid ar- 
rows. Each electron (hole) stripe in one layer faces a hole 
(electron) stripe in the other. The average filling factor of 
the highest occupied Landau level is v while that in the top 
layer is 1 — v, giving a total filling factor vt = 1. In the 
convention used here, a rung pair consists of the edges of 
an electron stripe in the top layer and an hole stripe in the 
bottom layer. The right-moving and left-moving chiral quasi- 
particle branches of each element of a rung pair are local- 
ized in opposite layers and denoted by A = 1,2. The mo- 
mentum-conserving back-scattering interactions not present 
in the Luttinger liquid model, discussed in this section, in- 
clude both interlayer and intralayer processes which have 
different behavior. The figure illustrates interlayer and in- 
tralayer back-scattering processes with the smallest possible 
momentum transfer. 

forward scattering processes only. Backscattering pro- 
cesses involve large momentum transfer and their bare 
matrix elements will be smaller in magnitude (see below) . 
We treat their effects perturbatively, using a renormal- 
ization group approach to account for the infrared di- 
vergences ubiquitous in quasi one-dimensional electron 
systems. The smectic state is stable only if the back- 
scattering interactions are irrelevant. 

The classical quadratic Hamiltonian which describes 
the energetics of small stripe edge fluctuations has the 
following general form 19 : 

^ = ^i dx i dx ' £ E E 

J J j,k=-oc A,ju=l,2 a,l3=R,L 

u$ a (x)K^(x-x',j-k)u% (x')] , (1) 

where the indices j and k label rung pairs, A and n the 
two different rungs within a rung pair, a and (3 the right 
or left moving chiral edges within a single rung, and x, 
x 1 positions along the stripes. In this equation, u^ a is 
the transverse displacement of the edge (j, A, a) from its 
classical ground state location. 

In Eq. (1), the linear charge density associated with 
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an edge displacement is p$ a (x) 
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the two-dimensional electron density inside the stripes 
n = l/2nl 2 and a = — , +, or L,R, for the left, right, 
moving fermions respectively. It follows from symmetry 
considerations that the elastic kernel satisfies the follow- 
ing equalities, 



(x), where n is p$ R (x) = : (x)R$(x) := R^(x)R$(x) - (rV(x)R$(x)) , (7) 



K%U) = K ZU) = (J) - . (2) 

which allow the Hamiltonian to be rewritten into sums 
and differences of the positions of left and right-going 
branch edges: 



T~to 



1 

8F 



dx 



^' E E 

j,k—-~oQ A, £i=l, 2 



(3) 



[u£*(x) + ul L {x)] K?{x - x',j - k) + < L (x')] + 

[ul R {x) - u$, L (x)] K*»{x - x',j ~ k) [u^ R (x>) - ul L {x')] 
where 



K^(x,j)=£ 2 ^K^(x,j) 

a/3 



a/3 



(4) 



(5) 



We regard the first term in this Hamiltonian as the con- 
tribution from fluctuations in the position of the rungs 
while the second term is the contribution from fluctua- 
tions of their charge densities. (The total filling factor 
varies locally when left and right going branch edges do 
not move together.) These two terms are analogous re- 
spectively to current and charge terms in the effective 
Hamiltonian of a conventional one-dimensional electron 
gas. The calculations we perform here will require only 
the long-wavelength limits of the x — x' dependence of 
elastic kernel in this Hamiltonian, which we estimate us- 
ing a weak-coupling approximation that we discuss be- 
low. 



B. Bosonization 

This Hamiltonian is quantized by recognizing that 
charge and position fluctuations result from particle-hole 
excitations at the edges of chiral quasiparticle branches, 
just as in an ordinary one-dimensional electron system. 
The real spin is frozen due to the presence of strong per- 
pendicular magnetic field and as a result we bosonize 
according to spinless bosonization scheme 21 . It follows 
from standard arguments that 

pl a {x),PkA x '>\ = ^ S ^ 6 ^ 3 .kd x S(x - x') . (6) 
In terms of Fermion creation and annihilation operators 



p x jL {x) = : L)\x)L${x) := Lf{x)L^{x) - {l)\x)L*{x)) , (8) 

with A G {1,2} denoting the rung in rung pair j and R, 
L, labeling right and left movers at the stripe edges. 

The low energy Hamiltonian is more conveniently de- 
scribed in terms of boson fields. The right and left mov- 
ing fermionic fields on the left stripe edge of rung pair j 
are given by 



while those on the right are given by, 



(9) 



(10) 



The above equations hold similarly for the left movers 
with the only change R — > L. Here b = a /I 2 is the 
width in k space of a rung and kp = av /2l 2 is the Fermi 
wavevector for the bottom layer stripes. The right and 
left slow fields can be expressed in terms of boson fields 
as in conventional one-dimensional electron systems: 



R^{x) 



2tt 



L*(x) 



1 



2tt 



(11) 



where <j>^ R {x) and (j>^ L {x) are the chiral components of 

the bosonic field $$(x) = [<f>j tR (x) + 4>j X ( x )] A 

In terms of the bosonic fields the chiral currents take 
the following form 



a 
2^ 



d x ^ a (x) 



(12) 



Introducing the dual field Qj(x) = [<pj, R (x) - <pj. L ( x )] A 
the position and charge variables U^, U'j, of the two 
edge system can be expressed as Uj — —l 2 d x $j and 

U'j = -l 2 d x G$. This shows that the field $ is related to 
the position fluctuations of the two edge system whereas 
is related with their charge density fluctuations. 

The field <& and the dual field 6 satisfy the following 
commutation relation 

[Qj(x),d x >^(x')} = -iwS j<k S x<ft 8(x-x') . (13) 

(The fields Q/^/n and —d x ^j\pK are canonical conju- 
gates). In terms of these new fields the Hamiltonian takes 
the following form 

Ho = \ fdxdx' J2[d x $j(x)K^(x - x',j - k)d x ,$£(x')+ 



3,k 



+d x Q x Ax)K^{x - x',j - k)d x ,Q£(x>) 



(14) 



whereas the Fourier transform of the corresponding ac- 
tion 
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3 A 



reads 



So 



^U*'(q lW )^(q)^(q, W ) 



A,p 



+^e A *(q^)^(q)e^(q,c) 



Here we have employed the shorthand notation, 



-A 

q,u> J —A 



dq x f" dq y a 
2tt _k 2tt 
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(16) 



(17) 



with A ~ l/l a high momentum cutoff, and have adopted 
the following Fourier transform conventions: 



Fj(x,t)= f e i ( q " x+q * a >- UT )F(q,L>) , 
F(q,w) = /* dxdr^e-^^^^'-^F^x,^ 



(18) 



(19) 



The QFT applies for distances larger than I and as a 
result the q x integration has to be cut off by ±2n/l. In 
Eq. (16) the kernel matrices K&(q) and K@(q) are the 
Fourier transforms of 



the consequences of weak quantum fluctuations when the 
ground state is similar to the mean-ficld-thcory stripe 
state. In this spirit, we use weak-coupling expressions 
for the interaction parameters of the model, replacing 
scattering amplitudes by the bare values for scattering 
of the Hartrcc-Fock theory quasiparticles. If the true 
ground state were a smectic, the values of these param- 
eters would be renormalized somewhat by higher order 
scattering processes. We expect that it will prove difficult 
to systematically improve on the estimates given here for 
the quantum Hall bilayer case because of the absence of a 
one-body kinetic energy term in the relevant microscopic 
Hamiltonian that would enable a systematic perturbative 
expansion. We emphasize that a quantitative theory of 
the forward scattering amplitudes that has a sound mi- 
croscopic foundation is necessary in order to decide on 
the relevance of the back-scattering interactions we have 
neglected so far and the character of the true ground 
state. As emphasized by the work of Fradkin, Kivel- 
son and co-workers 27 any conclusion is possible if the 
forward scattering interactions are allowed to vary arbi- 
trarily. The perturbative renormalization group scaling 
dimensions that we evaluate below are dependent only 
on the elastic constants at q x = 0, i.e. for strait stripe 
edges. The weak-fluctuation Hamiltonian may be evalu- 
ated in this limit by calculating the expectation value of 
the microscopic Hamiltonian in the Hartree-Fock theory 
ground state, which in this limit is a single Slater deter- 
minant with straight stripe edges displaced from those in 
the Hartree-Fock theory stripe ground state. By evaluat- 
ing the expectation value of the Microscopic Hamiltonian 
in a state with arbitrary stripe edge locations we find that 
for j ^ 0: 



K$»(x,j) = 2f[K^ R (x,j)+K^ L (x,j)] , (20) J dxK ^ /@{xJ) = 



K&{xJ)=Ky(xJ) = 2(?[K&{xJ)-K&{x,jj\ . (21) 

Integration over the <& fields in (16) yields an effective 
action in terms of the 8 fields alone 
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V(ja) T W(ja) W T V(ja - au) 
WTV(ja + av) V(ja) T W(ja) 



(23) 



e A *(q,u;) 



• (22) 



In the off diagonal elements the argument of W is the 
same as that of V, (ja±av) 1 respectively. The V and W 
contributions are proportional to two-particle intralayer 
and interlayer interaction matrix elements respectively, 
and are given by 



The corresponding 5$ action, obtained by integrating 
out the 6 fields, differs only through the interchange of 
<!> and 9, and K$, Kq. 



V(Y) 



■J 



V s N (q)e 



-iqY 



e 2- 



dq 



V S N (y/YW+f) , (24) 



C. Microscopic Theory of Long- Wavelength 
Interaction Parameters 

The objective of this coupled Luttinger liquid model 
for stripe states in quantum Hall bilayers is to address 



W(X) - J ^e-^ 2 V»(q)e-"> Y . (25) 

Note that the intralayer interaction contributions have 
competing direct and exchange contributions that cancel 



•5 



for Y = whereas the interlayer interaction has only 
a direct contribution. In the following we shall assume 
infinitely narrow quantum wells in both layers so that the 
interaction potentials occurring in the above equations 
read 



L 



N 



V'wd(?) , (26) 



where Ln(x) is the Laguerre polynomial form factor for 
electrons in the iV-th excited Landau level, and Vgj D is 
the Fourier transform of the Coulomb interaction within 
and between the layers, Vg(q) = e 2 27r/|g|, V^(q) = 
e 2 (27r/|g|) exp(— \qd\) with d being the layer separation. 
The long-ranged nature of the Coulomb interaction leads 
to logarithmic divergences in V and W which we reg- 
ularize by adding a term — (e 2 2-7r/|g|) exp(— 2|g|d ga t e ) to 
Vg, D with d gatc d. This regularization can be roughly 
thought of as introducing a metallic screening plane at 
distance d ga t leading to image charges that screen in- 
teractions between electrons in the bilayer system. Al- 
though V and W diverge for d gatc — > oo, it is possible to 
show 28 that K$, and Kq remain finite. In the following 
we choose a large but finite value for d ga t e for numerical 
convenience. 

The above form of the smectic energy kernel it$/e(j) 
applies for j ^ 0. For j = the components K R l(0) = 
K LR (0) and of K^ R (0) = K^ L (0) for A ^ ^ are given by 
the same expressions. The quantities K^ R (Q) — Kj} L (0) 
[= ^1^(0) = K 22 L (0)] have additional contributions 
that originate from the wavevector dependence of the 
Hartree-Fock self-energy at a given stripe edge and cap- 
ture the key property that the energy of the smectic 
must be invariant under rigid translations of all stripes, 
u ja( x ) l— * u ja( x ) + constant 19,29 . We find that 



K R \{q y ) 



k r 2 M 



D. The Balanced Bilayer Limit 

The special case of half filling in each layer has addi- 
tional symmetry that is most conveniently exploited by 
taking a slightly different approach. In this case the elec- 
tron stripes in both layers have the same width and the 
system therefore has effective periodicity of a/2. To be 
more precise, the problem can be formulated as a one- 
dimensional lattice of equidistant double edges placed a 
distance a/2 apart, noting carefully that right and left 
goers interchange their layer labels on alternate edges. 
To describe this system instead of using coupling ma- 
trices K^(x,j) and Kq L (x, j), as in the unbiased case, 
one can simply use the coupling constants Kq>(x,j) and 
Ke(x,j) 



K*,e{x,j) = 2l 2 [K RR (x,j)±K RL (x,j)} 



(30) 



with the value for K RR (0) reflecting translation invari- 
ance 



K RR (0) 



1 



+ 00 
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E (-^[V(j)-W(j)} . (31) 



In momentum space K§ and Kq> have the following form 



"<*-T>-"<T> 



■W(q, 



a a 



K e (q y ) 
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V(q y )-V(-) + W(q y ) + W(-) 



(32) 



• (33) 



We have used this simpler and partially independent for- 
mulation of the v = 1/2 limit, to test our results for the 
general case. 



<iy =0 



Note that these properties imply that det[Kg,(q y = 0)] = 
0. When these long wavelength approximations are 
employed, the Fourier transforms of K^, & in (23) de- 
pend only on q y but not on q x . From the relation 
K<s,/e(x,-j) = Kl /@ (x,j) it follows that K< s ,/ @ (q y ) is 
hcrmitian. Under particle-hole transformation, v i— » 
1 — v, the diagonal elements of these matrices remain 
unchanged while the off-diagonal elements transform as 

Ki%(ly) - e-^ a Kl) @ (q y ) , (28) 
Kl) e (Qy) ~ e +i ^Ki 2 (q y ) . (29) 



III. SMECTIC STATE PROPERTIES 

A peculiar property of quantum Hall stripe states is 
that the microscopic scale of back-scattering interactions 
is weak. For this reason observable properties may be 
those of smectic states over a wide interval of tempera- 
ture, even when back-scattering interactions are relevant 
at the smectic fixed point. In this section we discuss 
some characteristic properties of quantum Hall bilayer 
stripe states. 

A. Collective modes and Thermodynamic Properties 

The coupled Luttinger liquid model for bilayer quan- 
tum Hall stripe states gives rise to two collective modes 
with dispersions that can be determined by evaluating 
zeros of the determinant of the 2x2 matrix that defines 
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FIG. 2. Collective modes of the bilayer QH smectic phase. 
These results shown here were evaluated using /q (x, j) val- 
ues from Eqs. (23-25) with N = 2, v » 1/2 and a = 5.81. 
Notice that cu+(q x ,q y ) always disperses linearly for small q x , 
whereas u>-(q x ,q y ) disperses sub-linearly (as q x ) when q y = 0. 
The q y = behavior of the lower collective mode is sensitive 
to the q x dependence of the interaction coefficients which we 
do not evaluate microscopically. This illustration was con- 
structed by adding a small q x contribution to the interaction 
coefficients. 

the real time quadratic action at each q and u>. Writ- 
ing this matrix (with indices suppressed) as K^ 1 * 
[u 2 + 7T 2 <^X$ife]/?T 2 , it follows that the squares of the 
quadratic boson collective mode energies are 
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Both modes have energies that are proportional to q x . 
The velocity of the w_(q) mode vanishes for q y — » 0. 
In fact, when the q x dependence of K$ and Kq is 
dropped, as in most of our calculations u>- (q x ,q y = 0) 
vanishes identically; when the q x dependences are re- 
stored w 2 (q) ~ Qxily + Qx) at small wavectors and 
w - (Qx iQy = 0) oc q^ . For gate-screened Coulomb interac- 
tions, the x-direction w_ mode velocity is proportional to 
\q y \ in the small q x and q y limit. In the independent layer 
limit, the two modes become degenerate and we recover 
the isolated layer results obtained previously 19 . 

In the case of balanced bilayers the alternate formula- 
tion mentioned above is more convenient. The collective 
modes for this limit may be expressed as 



wi,2(q) = nq x \J K^(q)K e (q) 



(35) 



where K$(q) and ife(q) are given by Eqs. (32,33). The 
two collective modes of the general formulation applied to 
the v = 1/2 case correspond to two different wavevectors 
of this dispersion relation. 

The collective modes of the bilayer QH smectic phase 
are shown in Fig. 2. The right panel shows ui+(q x ,q y ) 
which disperses linearly in small q x for arbitrary q y . In 
contrast, cj_ (q x , q y ) disperses linearly at small q x for q y ^ 
0, but for q y = it is sub-linear: u)-(q x ,q y = 0) ~ q x . 

The Thermodynamic properties of the smectic phase 
of the bilayer system are those of a non-interacting bo- 
son system and are readily evaluated given the collective 



mode energies. For example, for v = 1/2 in each layer, 
using the simpler alternative formulation, we have one 
collective mode at each wavevector. The internal energy 
density is 
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where w(q) is given by Eq. (35). At low temperatures 
only the long-wavelength behavior matters and we obtain 
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where fe^To = y Kg(0)K$(0)/al. The specific heat of 
this system varies as T\u(Tq/T) at small T, vanishing 
less quickly than that of a non-interacting Fermion sys- 
tem because of the soft collective modes that result from 
the translational invariance of the stripe state. This low 
temperature behavior reflects the form of the dispersion 
relation at small q, w_(q) ~ Q x Qy\ only the prefactor of 
this enhanced specific heat changes in the unbalanced bi- 
layer case. These results for the specific heat are similar 
to the ones obtained in previous works by Barci et al. 24: , 
and Lopatnikova et al. 23 for the case of a single layer. 
There is no qualitative difference between the thermody- 
namic properties of single layer and bilayer stripe states. 
For unbalanced bilayers the specific heat at low temper- 
atures is dominated by the softer of the two collective 
modes, whose long wavelength dispersion is given by 



w-(q) ~ q x q y \ 



Tr[K«(p)K e (0)] 



det[#e(0)][det (#*)]"(()) 
It follows that the internal energy is given by 



(38) 



Tr [K*(p)Kg(p)] 



T 2 C(2)ln(^) , (39) 



Tr 3 V det(^(0))[det(^)]"(0)" 



and the specific heat will vary again as T\u(Tq/T). The 
To in Eq. (39) and Eq. (37) are given by corresponding 
expressions. 



B. Boson and Fermion correlation functions at the 
smectic fixed point 

In this section we discuss the static and dynamic cor- 
relation functions of the right and left moving fermions 
of the stripe edges and the Boson correlation functions 
in terms of which they are evaluated. The right and left 
moving fermion fields are expressed in terms of the $ and 
boson fields as follows: 
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where the $ held is related to position fluctuations of the 
two edge system, while the O field is related to its charge 
density fluctuations. We observe in the following that 
the charge and position fluctuations of the edges have a 
dramatically different effect on the correlation functions 
of the right and left movers. The single-particle Green's 
function for the right movers is given by 



i? At (x,r)i? A (0,0) 



(42) 
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We first evaluate the $ and Q field correlation functions 
C$ x and Cq X , where 



This result is generalized here by the average over q y 
and by the particular way in which the matrix nature of 
the Kq and K$ expressions enter the matrix elements 
above. The result for the field correlation functions 
differs only by the interchange of the Kq and K$ matri- 
ces. At hrst sight it appears that the position fluctuation 
factor in the right mover correlation function decays al- 
gebraically along the stripes. However, the power which 
characterizes this decay, d$. Xl is given by the integral 
over q y of Eq (46), which diverges logarithmically be- 
cause K- cx i>_ vanishes | — > 0; the same 
soft position fluctuations that led above to an enhanced 
specific heat, lead here to fermion correlation functions 
that decay faster than any power low but slower than 
an exponential. This observation generalizes to bilayers, 
a property of single layer stripe states noted by Lopat- 
nikova et aP 3 and Barci et aP 4 . 

Cq(x, 0,0), which specifies the charge fluctuation fac- 
tor in the fermion correlation functions, is given by 

C @ x (x,0,0)* (47) 



C xx (x,0,0) = ([$ x (x,t)-$ x (0,t)]- 



(43) 



and similarly for Cq X . In Eq. (43) and in the following 
the arguments of C*^ A e are (x — x',j — j',r — r'). From 
Eq. (22) we have for the $ A field 



C^ A (x,0,0) = 2 



where 
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Me(q,w) can be obtained by interchanging K<s>, Kq. 
The integral over ui is readily valued by decomposing 
M^ 1 as a sum over eigenmode contributions, writing it 
in the form ^ ± C+ 1 / [u> 2 + w±(q)] . ft follows after some 
algebra that correlation function can be expressed in the 
form 
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The charge fluctuation factor in the fermion correlation 
functions has a conventional algebraic decay with finite 
power d@ tX . The faster than algebraic decay of the 
fermion one-particle Green's function implies that the 
singularity in Landau gauge occupation numbers, a step 
( z tt)iction of unit magnitude in Hartree-Fock theory, is ex- 
ceedingly weak. 

The correlation function of the $ field along directions 
perpendicular to the stripes is given by 

'■ V dq x r' a dq y 

Qx J-Tr/a „ 



AAA 



(0,»,0) 



1/L 
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K+ + K- 



K xx 



K+K. 



(48) 



where C is a hnite constant that can be found numer- 
ically. It follows that the corresponding factor in the 
one-particle Green's function has a faster than algebraic 
decay in the thermodynamic limit. The factor associated 
with charge fluctuations has a similar dependence and is 
is given by 

T/ ° d Qy , 

— cos (%y) 



ono,2/,o) 



f A dfe r 

Jl/L 1x J- 



■/a 



for large x, where K± = v±(q)/ir. In the limit of bal- 
anced filling fraction for which Kq, and Kq are scalars 
the integrand which is averaged over q y above reduces 

to \J Kq I K§ , the familiar result for $-ficld correlation 
functions in a standard one-dimensional electron system. 



det(K^) 



K+ + K. 
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] + K+K. 
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where Cq is a function of the ratio (a/\y\) and is finite 
as \y\ — > +00. 



Similar steps can be taken for the general case of un- 
balanced bilayers. After the w integration, for (51) we 
obtain 



C. Tunneling density of States C<i>(0, 0, r) 

These results for boson correlation functions may be 
assembled to evaluate the imaginary time dependence of <{ | K 4 

the local fermion Matsubara Green's function and, by 
Laplace transforming this, the density of states for tun- 
neling into the bilayer system, a quantity that is in prin- 
ciple measurable. The single-particle Matsubara Green's 
function is given by 

(7(0, 0,r) = (i? A ](x,r)i? J A (x,0)) «cxp{-(l/2)C(r)} , (50) 
where (7(0,0, r) = <7$(0,0,t) + C e (0,0,r) and 
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We first discuss the balanced bilayer case for which the 
K^ 1 and Kq 1 matrices become simple numbers and the 
integral is simpler to treat analytically. In this case 
K<s,{q y ), K@(q y ), are given by (32-33) and 



Since at small q y , K_ ~ j^l, the most important contri- 
butions to the integral will come from the terms contain- 
\ ing 1/K- and from the exponential factors containing the 
argument u-(q y ). Keeping only these terms we obtain 
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We can understand the content of this integral by means 
of the following analysis. The integral can be separated 
into the sum of two terms, contributions from the region 
where q x and q y are small and the exponential can be 
approximated by the first few terms of the Taylor expan- 
sion and contributions from larger q x and q y where the 
exponential can be disregarded. The leading contribu- 
tion to the integral comes form the lower boundary of 
the second region where q y w l/(q x r). We focus on the 
case of large r for which our low energy theory applies. 
In this limit we can approximate KQ(q y ) w K@(0) and 
K<s,(q y ) w K<p(0)q y . In this limit Eq. (52) becomes 



demonstrating that the form of the Matsubara Green's 
function does not change qualitatively at unbalanced fill- 
ing factors. The tunneling density of states is the inverse 
Laplace transform of (7(0, 0, t), 



r+oo 
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where a and flo can be identified from Eq. (54) for bal- 
anced bilayers and from Eq. (56) in the unbalanced case. 
The inverse Laplace transform of such a function is 
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where Tfc(l) is the fcth derivative of l/r(r) at r = 1 and 
Hk arc the Hermite polynomials. In the asymptotic case 
of small E (low energies) we have 




*e(0) ln 2 



al 

The Matsubara Green's function factor contributed by 
field correlations has a weaker r dependence which we 
can neglect for the present qualitative analysis. 
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This shows that the density of states vanishes at the fermi 
energy stronger than any power of E and that in this 
sense the stripe edge physics in bilayer quantum Hall 
systems is not that of a usual system of weakly coupled 
Luttinger liquids. The above result generalizes to the 
case of bilayer systems, points that have been made about 
single-layer quantum Hall stripe states by Lopatnikova et 
aP 3 and Barci et aP 4 . 



IV. STABILITY OF THE QUANTUM HALL 
SMECTIC PHASE 

We now consider "backward" (interchannel) scattering 
interactions that do not conserve the number of electrons 
in each stripe edge. The most important conclusion of 
the following analysis is that back-scattering interactions 
are much more important for bilayer stripe states than 
for single-layer stripe states. We can classify the back- 
scattering interactions as either intralayer interactions 
that involve electrons only in one layer or interlayer in- 
teractions that involve electrons in both layers. Non-zero 
back-scattering two-particle matrix elements conserve to- 
tal momentum along the stripes, which means that the 
two Landau gauge guiding center jumps must sum to 
zero. In the case of quantum Hall stripe states, the 
microcopic matrix elements associated with these back- 
scattering processes tend to be small and these interac- 
tions will be important only if they produce strong in- 
frared divergences in perturbation theory. The strength 
of these divergences is characterized here by evaluating 
lowest order perturbative renormalization group scaling 
dimensions for these operators. As in Ref. 19 our renor- 
malization group (RG) scheme involves only x and r di- 
mensions and and treats the rung label as an internal 
index of the fields Qj. The philosophy underlying 
this procedure is discussed in the next section. 



A. Interlayer Tunneling 

We first address the tunneling of an electron from one 
stripe edge to the other in the same rung. The action con- 
tribution from this process has the following bosonized 
form: 



can be tested by considering the lowest order RG flow 
equation, 



->Tunn 



', [ dxdr [exp(i26^) + h.c] , 



(59) 



where the microscopic amplitude u is discussed below. 
We integrate out "fast" boson modes $ A , A in a shell, 
with A/6 < \q x \ < A and u),q y unrestricted, and then 
rescale q' x = bq x and uj' = bui leaving q y unchanged. With 
an appropriate rescaling of <f>, this RG transformation 
leaves the harmonic smectic action So [and the disper- 
sion relation (34) for -ftT$/e (23)], invariant. Stability of 
the smectic fixed point in the presence of back-scattering 



du 

= (2 - A T unn)M 



(60) 



with t = In b. As can be seen from this equation the 
tunneling operator will become relevant when its scal- 
ing dimension is less than two. When this term in the 
Hamiltonian is strong, the system is described at low en- 
ergies by a set of quantum sine-Gordon models coupled 
by gradient terms (22). In the domain < Axunn < 2 
the continuous symmetry, present in (22), and broken by 
tunneling, is lost in the low-energy fixed point action. 
This model has a discrete symmetry A — > 6 A + nn for 
any integer n and the QFT becomes massive 30 . The gap 
due to tunneling will lead to an integer quantum Hall 
effect at total filling factor v T — I. Using Eqs. (22), 
and (59) , we find the following expression for the scaling 
dimension 
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r /a dqy 

/ 27T 
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The integrand in the integral over q y is similar to that 
involved in the boson correlation field and, ignoring 
the matrix character of the coefficients th at appea r in the 
smectic fixed point Hamiltonian, is ~ W Kg / Kq . Since 
Kg vanishes for q y — ► 0, we can expect this quantity 
to be small. Indeed we find by evaluating this integral 
numerically that interlayer tunneling is always relevant. 



B. Coulomb Backscattering Interactions 

The tunneling amplitude in bilayer quantum Hall sys- 
tems can be made extremely small by making the bar- 
rier between quantum wells higher or wider and is often 
completely negligible in practice. Coulomb interactions, 
on the other hand, are always present and must always 
be considered. We consider interlayer and intralayer 
Coulomb back-scattering processes separately. In the 
strongest interlayer back-scattering process an electron is 
transferred from, say, a left-moving top layer stripe edge 
to a right-moving edge in the same rung pair of the same 
layer, while in the same rung pair of the bottom layer 
an electron is transferred in the opposite direction, as 
depicted in Fig. 1. The interlayer back-scattering opera- 
tors for processes involving neighboring rungs have large 
scaling dimensions and tend to be irrelevant. In addi- 
tion, the bare matrix elements for such a process will 
fall off rapidly in magnitude with increasing distance be- 
tween the rungs involved. The action for interlayer back- 
scattering interactions for electrons within the same rung 
pair reads 
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s J dxdr [exp(i2(6] + 6 2 )) + h.c] . Ae = T 7 " ^% 

3 J-Tv/a -T 



det K$ 



(1 - cos(q y a)) 



(62) 

[The other kind of process involving two neighboring rung 
pairs is related to the above one by a particle-hole trans- 
formation v i— > 1 — v\ . 

After an elementary calculation we obtain the following 
scaling dimension expression: 



2 det K* 




■ (63) 



This expression is similar to that which would be ob- 
tained for inter- wire back-scattering interactions in a sys- 
tems of two coupled quantum wires. This integral is 
similar to the one that appears in the tunneling oper- 
ator scaling dimension calculation, although it is easy to 
verify that forward scattering interactions between dif- 
ferent stripes play an essential role. As we discuss below, 
this operator is usually strongly relevant (Ai n t or — » 0), 
so that at low temperatures the phases 1 and 9 2 of 
neighboring two edge system are strongly anti corre- 
lated. The low energy nontopological (chargeless) exci- 
tations in this limit can be understood by approximating 
cos[(9} + 9 2 )] ss 1 - (6} + 6 2 ) 2 /2. When a term of this 
form is added to the quadratic Hamiltonian, the low- 
energy collective mode dispersion at long-wavelengths 
takes the form of a spatially anisotropic two-dimensional 
XY ferromagnet, uo 2 <~ Kq 2 + uq 2 . We discuss the signif- 
icance of this result at greater length below. 

Finally, in an intra-layer back-scattering process two 
electrons move in opposite directions between pairs of 
stripe edges in the same layer with the same separation 
(cf. Fig 1). Here we also concentrate on processes in- 
volving neighboring rungs only. Processes involving two 
rung pairs are again related to those involving three pairs 
by a particle-hole transformation. For the first case the 
action reads 

Sintra = -« / dxdr [exp(i($ 2 - $) + - $ 2 +1 )) 

• exp^-e 2 - e 1 + e) +1 + e 2 +1 )) + h.c] , (64) 

which leads to a scaling dimension Ai ntra = A$ + Ae 
with 
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K 11 



K 12 
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Note that the imaginary part of the integrand in 
Eqs. (63), (65-66) does not contribute to the integrals. 
Backscattering processes other than the ones discussed 
above have larger scaling dimensions and also involve 
larger momentum transfer and have therefore exponen- 
tially smaller bare matrix elements. We therefore shall 
concentrate on the processes discussed above. 

We now discuss our numerical results for scaling di- 
mensions of the operators that are not described by the 
quadratic boson theory. Fig. 3 shows the scaling dimen- 
sions of the back-scattering interactions in a balanced 
system (y = 1/2) in the second excited Landau level 
(N = 2) as a function of the layer separation djl. The 
stripe period chosen for these calculations, a = 5.8£, 
corresponds to the period at which the Hartrec-Fock 
energy of the stripe state in a isolated layer is mini- 
mized. 31 . Interestingly, single-electron tunneling is ir- 
relevant (Axunn > 2) for d/£ > 1.5, but is strongly 
relevant at smaller layer separations. Interlayer back- 
scattering is relevant (Ai ntcr < 2) at all layer separa- 
tions, more strongly so at smaller layer separations, while 
the scaling dimension of the intra layer back-scattering is 
smaller than 2 only for d/£ 2 and approaches a value 
of A intra ~ 1-84 for d/t > 1. This value for the limit of 
weak interactions between the layers recovers the single- 
layer result obtained earlier. 19 . The contributions A$ 
and Ae to Ai ntra , not shown in the figure, become equal 
in this case. 




4 6 8 10 

layer separation d/l 

FIG. 3. Scaling dimensions for tunneling and back scat- 
tering interactions as a function of layer separation d/l in a 
balanced bilayer system at bilayer total filling factor vt = 9, 
i.e. with a N — 2 valence Landau level. 
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FIG. 4. The left panel shows the N = 2, d/t = 5.8, in- 
tralayer back-scattering interaction scaling dimension. The 
right panel shows the scaling dimension of the Tunneling op- 
erator. 




in " r <U,=50.0 N= 2.0 a/i= 5.i 




FIG. 5. The left panel show the scaling dimensions for in- 
terlayer back-scattering across narrow and the right panel for 
interlayer back-scattering across wide rungs in unbalanced bi- 
layers. Note the difference in scale between left and right pan- 
els. The most relevant interlayer back-scattering interactions 
are those of narrow rungs. 

The dependence of scaling dimensions on the bilayer 
balance is illustrated in Fig. 4 and Fig. 5. We see that in- 
tralayer interactions become more relevant when their in- 
dividual filling factors move away from v = 0.5, as in the 
single-layer case, while the tunneling operator becomes 
less relevant. Interestingly the interlayer back-scattering 
interactions show different results depending on the dis- 
tance between the edges involved in the transition. For 
v 7^ 0.5 we have to distinguish between nearest neighbor 
interlayer and intralayer back-scattering processes that 
involve, according to the definition given in Fig. 1, only 
the smallest number of neighboring rung pairs (one and 
two, respectively) and those processes that involve for- 
mally two and three rung pairs, respectively. These two 
kinds of processes are related by particle-hole transforma- 
tion and therefore shown in different panels. Generally 
the scaling dimension increases with the distance between 
the edges. The data shows that one of these two back- 
scattering processes, related by a particle-hole transfor- 
mation, is always relevant and that the minimum scaling 
dimension decreasing with increasing bias between the 
layers. In summary, the most relevant residual interac- 
tions is interlayer back-scattering and they are increas- 
ingly important as the bilayer is unbalanced. 

We note that the scaling dimensions ATunn of the 
single-electron tunneling and Ai nt0 r of the interlayer 
back-scattering approach zero for d/t — > 0, i.e. these 
processes become strongly relevant. This is a natural re- 
sult since in this limit we recover the monolayer electron 
spin quantum Hall ferromagnct. This system is perfectly 
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J= 4 
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N=0 a=3.3l 
v=0.5 



N=2 a=5.8l 



N=1 a=4.4l 



...... A) 

A ; „, 



N=3 a=6.9l 



6 8 2 4 
layer separation d/l 



8 10 



FIG. 6. The scaling dimension of the back-scattering inter- 
actions in the various Landau level (N — 0, 1,2,3) in a bal- 
anced system as a function of the layer separation d/i. The 
stripe periods a are obtained from Hartree-Fock monolayer 
results given in Ref. 31 . 

isotropic in pseudospin space, and therefore processes 
like tunneling which acts essentially like a (pseudo- 
magnetic field are obviously very relevant. This in- 
creased relevance arises formally in our calculations 
through the property that the matrices K$(q y ) vanish 
in the limit d — > 0, so that the integrands in Eq. (61) 
and (63) are identically zero. Kq, vanishes because it 
is a measure of energy changes associated with charge 
transfer between layers at a particular stripe edge; when 
d — > only the total charge near each stripe edge influ- 
ences the energy functional. In Fig. 6 we show the depen- 
dence of these scaling dimensions on Landau level index 
N = 0, 1, 2, 3 with the stripe periods taken from Ref. 31 in 
each case. As this figure shows, our results for the scal- 
ing dimensions of back-scattering processes around the 
assumed stripe state depend only weakly on the Landau 
level index. We note that in the lowest and first excited 
Landau level (N — 0, 1) no conductance anisotropies are 
found experimentally in single layers, even though there 
is a stripe state in each of these Landau levels which is 
a local minimum of the Hartree-Fock energy functional. 
The true ground state in these instances is far from the 
stripe state, differencing in character even at microscopic 
length scales. The fact that our calculation does not ob- 
tain anomalous results in cases where we do not believe 
stripe states occur, emphasizes again that our approach 
can only address the properties of systems in which fluc- 
tuations around Hartree-Fock stripe states are weak. It 
cannot predict when stripe states occur. Future experi- 
mental activity will be necessary to identify with confi- 
dence when stripe states occur in bilayers. 
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C. Smectic interlayer phase coherent and Wigner 
Crystal states 

1. Smectic interlayer phase coherent state 

In this section we examine the effect of the interlayer 
back scattering interactions (when they are strongly rel- 
evant) on the low energy physics of the system and show 
that the phase coherence is marked by a nonvanishing 
value of an interlayer phase order parameter. In this 
phase electrons at each stripe edge are coherent superpo- 
sitions of the upper and lower layer states. 

The most relevant interlayer back-scattering operators 
are related by particle-hole symmetry and describe back- 
scattering across an electron stripe in the top layer and 
the corresponding hole stripe in the bottom layer or 
across a hole stripe in the top layer and the correspond- 
ing electron stripe in the bottom layer. In terms of the 
Luttingcr liquid fields we have defined, the sum of these 
two interactions takes the form: 



e 



Ointor = -Winter,!/ COs[2(©| 

-"inter.Cl-./) COs[2(9 2 + 9^ +1 )] 



(67) 



Expressions for the bare values of these coupling con- 
stants are given below. As shown on Figs. 5(a), 5(b), 
(left, right panel respectively), at small layer separations 
these operators are strongly relevant. At low tempera- 
tures the phases 9] and 9 2 , 9 2 and 9] +1 of neighboring 
two edges tend to be strongly anti correlated. The low 
energy excitations in this limit can be understood by ap- 
proximating cos[(6j + 9 2 )] w 1 - (9} + e2) 2 /2. When 
terms of this form are added to the action, it takes the 
following form: 



e A *(q, W )M*"e"(q,a;) . 

The new matrix Mo is given by 

u il +u i2 un + u i2 e~ ig v a 



(68) 



Me = M + 2 



Mil + u i2 e iq y a Mii+«i 2 



(69) 



where M@ is the matrix of the system at the smectic 
fixed point and is given by Eq. (22) or by Eq. (45) (inter- 
changing K<p, K@) and uu, u- l2l is the short notation for 
Winter^, Mi ntc r,(i-!/) > respectively. The effects of the inter- 
layer back-scattering interactions, included on the new 
matrix of Eq. (69) (which we denote by N@), shift the 
poles of the boson propagators. The low energy collective 
modes now are given by 
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(70) 
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FIG. 7. The collective modes of the bilayer QH smec- 
tic interlayer coherent phase, for the case of K&/q(x, j) of 
Eqs. (23-25) and N = 2, v « 1/2 and a = 5.8/, are shown. 
u)+(q x ,q y ), the right panel, is of the form of a spatially 
anisotropic two-dimensional ferromagnet, cj+(q) ~ Kq x +uq y . 
The u}-(q x ,q y ) collective mode vanishes only for q y — > and 
q x — > when non-local contributions to the interaction coef- 
ficients are accounted for. 

A = q 2 x Tr(K$K®) + 2N£K 22 + 23? {K 12 N£ } , (71) 

D = det(tf*) {q x det(Ke) + 2N 22 K^ + det(A e ) (72) 

-2&{K%Ng}} . 

These low energy collective modes are shown in Fig. 7. 

As before, the case of balanced filling fraction can be 
described in a more transparent way using the scalar cou- 
plings Kq of Eq. (32), (33). For this case the low- 
energy collective modes have the form 



<4(q) 



tt 2 #*(q) q 2 x K e (q) + 2u- a (1 + cos(q y a/2)) 



(73) 



In this formulation one of the gapless modes is located 
at the edge of the Brillouin zone, which is now doubled, 
at q y = 2tt /a. In the extended Brillouin zone we use for 
balanced bilayers the w_ (q) softmode appears for q y — ► 0, 
whereas the w+(q) softmode appears as q y — > 27r/a. The 
two modes have the following behaviors: 



wi(q) « n 2 Kl(0)ql K B {Q)q 2 x + 4un 



( 



^(q)«7r 2 ^*(^) 



Ke(-)q x + — 



27r\ 2 ' 

q y--) 



( 



Similar results can be obtained using the matrix formu- 
lation for general v and become equivalent for v = 1/2. 
There is no qualitative change in the collective mode 
structure when v ^ 1/2. 

The interlayer phase coherent smectic state is charac- 
terized by a finite value of the following order parameter, 



*(r) - (4(r)Vs(r) 



1 

27 



-2i0(r 



1 

2^ 



i(4e 2 (r)) 



where V'b are fermion creation and annihilation for 
the top and bottom layers respectively. We now show 



where 
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that (0 2 (r)) is finite. We discuss only the case of bal- 
anced bilayers, using the alternative formulation which 
is more transparent. We find that 

= 



/ 



dto 



Jf*(q) 



2ir 

a/2 



A"*(q) g2^ e ( q ) + 2uil [1 + cos(^)] 



— / 

167T J_ 



2-K/a 

dq yy — 

2?r/a V K e (q y ) 



In 



2^ e (g»)(l/0 
un(l + cos((7 y a/2)) 



2. Coherent smectic state specific heat 

The internal energy of the bilayer smectic phase coher- 
ent (SPC) state will be dominated at low energies by the 
contribution from the w_(q) mode. The leading contri- 
bution to the integral for the internal energy comes from 
the region of small q. The q x integral now has a natural 



infrared cut-off however at y Uu/K@(0). It follows that 
the internal energy is given for small u by 



(77) 



U: 



2a 



In (77) we have introduced an upper short distance cut- 
off l/£ for the q x integration. Interlayer back scattering 
interactions have cut-off the infrared divergence of the q x 
integration, making the integral finite and establishing 
particle-hole pair condensation. In this state the U(l) 
symmetry associated with conservation of total charge 
difference Nt — Nb between top and bottom layers is 
broken. 

We conclude on the basis of this analysis that interlayer 
back-scattering will drive the Hartree-Fock bilayer smec- 
tic state to a state which has both broken translational 
and oricntational symmetry and spontaneous interlayer 
phase coherence along the edges. We expect this state to 
exhibit giant interlayer tunneling conductance anomalies 
at low-bias voltages, similar to those that have been seen 
in the N = Landau level in bilayers. Although these 
states have a charge gap that we discuss below and should 
exhibit the quantum Hall effect, we expect that they will 
exhibit strongly anisotropic dissipative transport at finite 
temperatures. Their two gapless collective modes arise 
because they have broken translational and oricntational 
symmetry and spontaneous interlayer phase coherence. 
We also note that the quantum character of these bilayer 
smectic states is quite distinct from the quantum smec- 
tics discussed previously for the single-layer case. For 
instance the long-wavelength behavior of the quantized 
collective mode uj-(q x ,q y ) changes from being propor- 
tional to |<7a;9y| to being proportional to \q y \ only when 
spontaneous inter-layer phase coherence is present; lock- 
ing the phase difference between different layers qualita- 
tively increases the cost of independent position fluctu- 
ations. The long-wavelength behavior of the Lu+(q x ,q y ) 
collective mode is that of an anisotropic superfluid. As in 
the case of uniform states, spontaneous interlayer phase 
coherence is equivalent to electron-hole pair superfluid- 
ity, but the broken orientational symmetry of the smectic 
state causes this superfluid to have orientation dependent 
stiffness. 

It seems quite possible that the order parameter that 
characterizes the broken orientational and translational 
symmetry of these states will be driven to zero when in- 
terlayer interactions are sufficiently strong. Indeed this is 
suggested 15 ' 16 by mean-field calculations. We are unable 
to estimate where this transition takes place using the 
methods of this paper. 



7^1^(0)^(0) 



T 2 C(2) In 



K e (0) 1 
4iiii al 



(78) 



and the specific heat will now be linearly dependent on 
T. The specific heat anomaly noted previously for the 
bilayer smectic is suppressed when inter-layer coherence 
is established, even though broken translational and ori- 
entational symmetry are still present. 



3. Wigner crystal state 
Intralayer back-scattering interactions take the form 

Ointra = -u[cxp(l(2k F X + $ 2 - $} + - $ 2 +1 )) 

• exp( i (-e 2 -e 1 + e 1 +1 + e 2 +1 )) + h.c] , (79) 

where the oscillatory dependence on coordinate along the 
edge which we have exhibited explicitly follows from our 
earlier field operator definitions. This interaction domi- 
nates only at quite large layer separations. When it does 
it drives the system to a state which has periodicity along 
the stripe edges as well as across the stripes. Since, the 
wavelength along the stripe is 4nl 2 /a, and since the pe- 
riodicity along the direction perpendicular to the stripes 
is a, this state will contain one electron per layer per two 
dimensional unit cell. We therefore identify this state as 
a bilayer Wigner crystal (WC) state. 

D. Gap Estimates for Bilayer Stripe States 

The most important conclusion from the above calcu- 
lations is that interlayer Coulomb back-scattering inter- 
actions are always relevant in bilayer stripe states. The 
gapless bilayer stripe state can never be the true ground 
state. Since the bare matrix elements associated with 
these interactions are often quite small, however, they 
will often be important only at low temperatures. As we 
explain below, we believe that Coulomb interactions will 
most often drive the system either to an isotropic coher- 
ent state or to a smectic coherent state. Both states will 
have a charge gap and an integer quantum Hall effect. In 
this section we estimate the size of this gap and hence the 
temperature above which we expect the phenomenology 
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of these states to cross over from quantum Hall behavior 
to stripe state behavior. 

Our estimates are built on bare matrix elements whose 
evaluation we discuss below and on the scaling dimen- 
sion calculations discussed above. Given dimensionless 
inter- and intralayer back-scattering interactions u; n tor 
and Uintra, we can estimate the gap by integrating the 
RG flow equations to obtain 

E e g / a (u e/a )=b- 1 E^ a (b 2 - A ^u e/a ) , (80) 

where the super- and subscripts e/a are used for inter- 
and intralayer interactions respectively. When the inter- 
actions become of order 1 on the renormalized energy 
scale (b 2 ~ A u = 1), the energy gap should be roughly 
equal to the renormalized characteristic Coulomb energy 
E c , giving 

E g (u) = (U/E c ) 1 ^- A ^E c , (81) 

where U = uE c is the microscopic high-energy-scale 
back-scattering interaction strength. The v dependence 
of the gap enters through U, and through the scaling di- 
mensions. Both effects conspire to strongly reduce the 
gap magnitude near half filling. Taking E c = Q.3e 2 /l, 
approximately the maximum correlation energy per elec- 
tron in a partially filled Landau level, the resulting gaps 
for TV = 2 and d gat0 = 50Z, are shown as a function 
of filling fraction and distance between layers in Fig. 8 
and Fig. 10. We notice that the gap resulting from the 
intralayer back-scattering interaction is very small near 
half filling, dropping below the range accessible to dilu- 
tion fridges over most of the filling factor shown in this 
figure. On the other hand the gap resulting from the in- 
terlayer back-scattering interactions is not as small and 
remains reasonably large out to large values of the in- 
terlayer separation d. Recalling that this interaction is 
proportional to R\L\L\r 2i we see that when this inter- 
action is strong it favors interlayer phase coherence along 
each stripe edge and that when it is very strong it leads 
to condensation of the field Oj + <d 2 to a value indepen- 
dent of j. Since is by definition the phase difference 
between left and right-going fermion fields at the (j, A) 
stripe edge, and since the layers indices of right and left 
going fermions is opposite at A = 1 and A = 2 stripe 
edges, what is condensing when this interaction is strong 
is the phase difference between fermions in opposite lay- 
ers. In other words, the state that occurs in the strong 
interedge back-scattering limit has spontaneous interlayer 
phase coherence. States with interlayer phase coherence 
and stripe order can occur as local and even global min- 
ima of Hartree-Fock energy functionals. Coupled with 
the irrelevance of intralayer back scattering interactions 
at small layer separations in the bilayer case, our analy- 
sis suggests that they can be the ground states of bilayer 
quantum Hall systems in high Landau levels. 




FIG. 8. Estimated charge gap due to interlayer 
back-scattering interactions. These interactions are always 
relevant and lead, in the absence of interlayer tunneling, to 
states with spontaneous interlayer phase coherence. The en- 
ergy scale in this figure is ~ e 2 /e£ which is ~ fes50K for a 
typical higher Landau level experiment. The energies should 
be reduced to account for screening from inter-Landau level 
transitions that we have not included in our calculations. 



0.15 




0J ■ . ■ 1 ■ 1 ■ 1 ■ r- 

2 4 6 8 10 

d/l 



FIG. 9. Estimated charge gap due to interlayer 
back-scattering interactions, for balanced bilayers (y = 1/2 
in each layer), as a function of layer separation. This depen- 
dence is extracted from Fig. 8 and shown here for clarity. 

For intralayer back scattering, the bare backscattering 
interaction matrix element has both direct and exchange 
contributions, while interlayer back-scattering has only 
a direct contribution. An elementary calculation using 
Landau gauge basis states leads to the following explicit 
expressions that were used to obtain gap estimates: 
Same-Layer Direct 

hi 2 = Y 1 + Ql 2 , k 2 l 2 = Y 2 - Ql 2 , (82) 

k 3 l 2 = Y 1 , k 4 l 2 = Y 2 , (83) 
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FIG. 10. Estimated charge gap that would result from in- 
tralayer back scattering interactions in the bilayer case. When 
the scaling dimension is larger than two, the gap vanishes. 
Intralayer interactions are more important than interlayer in- 
teractions only at very large layer separations. The energy 
scale in this figure is e 2 /el. 

and 

(y 1 + qi 2 ,y 2 - Qz 2 |y|y 1; y 2 ) = ^ a2v2/212 ■ 

■ J dq x e-^ l2 / 2 e-^ a V s n (q x ,Q) , (84) 
Same-Layer Exchange 

k 2 l 2 =Yt + Ql 2 , k x l 2 = Y 2 - Ql 2 , (85) 

k 3 l 2 = Y 1 , k 4 l 2 = Y 2 , (86) 

(Y 2 - Ql 2 , Yi + Ql 2 \ V\Y U Y 2 ) = ^e- a ' 2 / 21 ' 2 ■ 

■ J dq x e-^ l2 / 2 e-^ a "V s n (q x ,a/l 2 ) , (87) 

Different-Layer Direct 

hi 2 = Yi + Ql 2 , k 2 l 2 = Yi , (88) 
k 3 l 2 = Y 1 , k i l 2 =Y 1 +Ql 2 , (89) 

(Y 2 ,Y 1 \V\Y 1 ,Y 2 ) = _L e -^ 2 /2i 2 . 

Ztt 

■ J dq x e-& 2 / 2 V£(q x ,Q) , (90) 

where the subscripts S and D refer to two-dimensional 
Fourier transforms of the Coulomb interactions between 
electrons in same and different layers. We see in Fig. 8 
and Fig. 10 that the importance of interlayer interactions 
diminishes rather slowly with layer separation, leading 
to sizable integer quantum Hall gaps out to large layer 
separations. 
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FIG. 11. Apparent phase diagram predicted for experimen- 
tal studies of high mobility bilayer systems at dilution fridge 
temperatures. The various phases in this illustration have 
qualitatively different transport properties. These calcula- 
tions are for stripe states in the N = 2 orbital Landau level 
(a — 5.8£), with weak remote-gate screening (d ga t c = 501). 
It is possible to explore the phase diagram experimentally in 
a single sample, since both the top layer filling factor v and 
the normalized interlayer separation d/£, are altered when 
the charge imbalance and total electron density are changed 
by using front and back gates in combination. For inter- 
layer spacing d less than approximately 1.5£ and any charge 
imbalance, we expect the bilayer to be in an isotropic in- 
terlayer phase coherent (IIPC) state which has a large gap, 
integer quantum Hall effect and isotropic transport proper- 
ties. Anisotropic states are expected only for more widely 
spaced layers, d > 1.5£. For strongly unbalanced layers (v far 
from 1/2) we expect anisotropic Wigner crystal (WC) states 
to appear because of intralayer backscattering interactions, 
just as they do in the single layer case. These states will 
exhibit a quantum Hall effect with an odd integer quantized 
Hall conductivity. Stripe (smectic metal) states (SS) tend to 
occur when each layer has a filling factor close to v = 1/2, 
but as in the single layer case these states are never the true 
ground states. Smectic metal states show anisotropic trans- 
port but do not show an integer quantum Hall effect. Inter- 
layer backscattering interactions always induce charge gaps 
but these are sometimes too small to be observable at a typical 
dilution fridge temperatures which we take to be 0.001e 2 /d. 
Regions with an estimated charge gap larger than this value 
are labeled as smectic phase coherent (SPC) state regions in 
the phase diagram. Smectic phase coherent states have an odd 
integer quantum Hall effect, and are expected to have trans- 
port properties which are much more anisotropic than those 
of the anisotropic Wigner crystal states. This state should 
also exhibit giant interlayer tunneling conductance anomalies 
at low bias voltages. 

Our results for the energy gaps are summarized in 
Fig. 1 1 by a schematic phase diagram intended to repre- 
sent predicted experimental findings in very high mobil- 
ity bilayer systems at dilution refrigerator temperatures. 
This phase diagram was constructed from a recipe spec- 
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ified below. Different regions of the phase diagram as 
a function of layer separation d/l and imbalance, char- 
acterized by are identified as exhibiting the behav- 
ior of one of the following phases. The bilayer smectic 
state is a state with no integer quantum Hall effect and 
anisotropic transport. The coherent bilayer smectic state 
will have an integer quantum Hall effect but will still have 
anisotropic transport at finite temperature. The bilayer 
Wigner crystal state will have an integer quantum Hall 
effect with an odd integer quantized Hall conductivity. 
We predict bilayer smectic state behavior when neither 
interlayer nor intralayer back-scattering interactions pro- 
duce a gap larger than 0.001e 2 /e£. We judge that a gap 
smaller than this size would not produce observable ef- 
fects in a typical dilution fridge experiments. Interlayer 
backscattering interactions are much more effective than 
intralayer interactions in producing gaps because they 
are strongly relevant. We predict bilayer Wigner crys- 
tal behavior when the intralayer back-scattering yields 
the largest gap and a gap that exceeds our minimum 
value. These states are expected only when the charge 
imbalance is large or the layer separation is quite large. 
We predict bilayer coherent smectic states when inter- 
layer back-scattering produces the largest gap, provided 
again that it exceeds our minimum value. Because the 
intralayer interactions are strongly relevant observable 
gaps are expected out to very large layer separations, 
an unexpected result of our analysis. The interval of 
charge imbalance where stripe (smectic metal) states are 
expected expands only modestly with layer separation, 
but is sensitive to the orbital index N of the Landau lev- 
els, since nodes in the orbital wavefunctions can cause the 
bare backscattering matrix element to vanish at partic- 
ular TV-dependent values of v. The details of boundaries 
separating stripe state and stripc-phasc-cohcrent regions 
of this phase diagram will be quite different for differ- 
ent values of N. As we have emphasized, our approach 
is reliable only when quantum fluctuations around the 
mean-field stripe state of Hartree Fock theory are weak. 
For small layer separations the charge gaps start to be- 
come comparable to the underlying microscopic energy 
scales. In this regime we expect that the ground state 
is actually an isotropic coherent bilayer state, but arc 
unable to provide a reliable quantitative estimate of the 
layer separation at which this transition occurs. 

V. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied double layer quantum 
Hall systems at odd integer total filling fractions. Mean 
field theory predicts that these systems can form striped 
ground states. This observation serves as the starting 
point for our work. The Hilbert space in which the low- 
energy excited states of mean-field bilayer stripe states 
reside may be mapped to those of an infinite set of cou- 
pled Luttinger liquids, one for each stripe, allowing us 



to borrow bosonization techniques from the literature on 
one-dimensional electron systems. Quantum fluctuations 
around the mean-field stripe state are conveniently de- 
scribed in terms of Bose quantum fields that can be in- 
terpreted as representing charge density and the position 
fluctuations along each stripe edge. The interactions that 
control quantum fluctuations in the electron ground state 
include both forward scattering terms which contribute 
to quadratic interactions in the Boson Hamiltonian and 
weak, but more complicated back-scattering terms. The 
coupled Luttinger liquid model obtained when the back- 
scattering interactions are neglected is not of the stan- 
dard form because both charge and position terms in 
the effective Hamiltonian have a matrix character and 
because the energy cost of fluctuations in which stripes 
move collectively is small when the stripes are not pinned. 
We find that the latter property leads to Fermion spatial 
correlations whose decay is faster than any power law, 
to a specific heat that vanishes less quickly than T for 
T — > 0, and to a tunneling density of states that vanishes 
faster than any power law for E — > 0. These properties 
of bilayer stripe states are similar to properties estab- 
lished previously for single layers by Lopatnikova et al. 
and Barci et al. There is no limit in which bilayer stripe 
quantum Hall states can be treated as a system of weakly 
coupled Luttinger liquids. 

We address the role played by intralayer and inter- 
layer back-scattering interactions by evaluating their per- 
turbative renormalization group scaling dimensions, fol- 
lowing an approach two of us have taken previously for 
the case of single- layer stripe states. 19 In the single-layer 
case we reached the conclusion that these interactions 
are always relevant, and that they likely drive the sys- 
tem to a Wigner crystal state with an energy gap. Esti- 
mates of the size of this gap based on bare back-scattering 
matrix elements and scaling dimensions gave extremely 
small values, however, consistent with the observation of 
stripe state phenomenology at temperature scales that 
could be reached experimentally Since other researchers 
have reached different conclusion about the relevance of 
back-scattering interactions in single layer systems, it is 
worthwhile in stating the conclusions we have reached 
in the present work to emphasize once again the philos- 
ophy that underpins our calculations and explain why 
we have considerable confidence in the conclusions we 
reached previously. 

Our identification of a low-energy Hilbert space in 
which it is possible to derive a simplified many-electron 
Hamiltonian is based on the experimental discovery of 
stripe states and on evidence from experiment that the 
true ground state is energetically very close to the mean- 
field theory ground state. In our view the most convinc- 
ing evidence in this regard is the ability 32 of Hartree- 
Fock theory to accurately predict the dependence of the 
stripe state orientation on in-plane field strength, quan- 
tum well width, and other microscopic parameters. In 
single-layer systems, quantum fluctuations arc impor- 
tant only at low-energies and long length scales. When 
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mean-field theory accurately describes the microscopic 
length scale physics, we can use the elementary exci- 
tations of the Hartree-Fock stripe state to identify the 
Hilbert space of low-energy excitations, and confidently 
use bare interaction matrix elements to estimate forward 
and back-scattering interaction parameters. The issue of 
quantum stability of smectic states in single-layer sys- 
tems has received interest partly because it is closely 
related to the possible existence 33 ' 27 ' 34 of freely sliding 
analogs of the Kosterlitz-Thouless phase in stacked two- 
dimensional XY models. Although it is certainly clear 27 
that interacting Luttingcr liquid fixed-point actions exist 
for which back-scattering interactions are irrelevant, that 
is not sufficient to decide on their relevance in the case 
of quantum Hall stripe states. Crudely speaking, irrel- 
evance in the case of repulsive interactions requires 19 ' 28 
that the forward scattering interaction strength decay in 
a strongly non-monotonic way with edge separation. For 
single-layer systems Fertig and collaborators 35 have esti- 
mated forward scattering amplitudes using an approach 
that goes beyond the weak coupling approximations we 
employ, doing so however in a partially ad hoc manner 
by fitting their model to collective modes evaluated in 
a time-dependent Hartree-Fock approximation. Their 
conclusion on the relevance of back-scattering interac- 
tions is opposite to ours. The source of the discrep- 
ancy may be traced to broken particle-hole symmetry 
in the half-filled Landau level Hartree-Fock approxima- 
tion Wigner crystal state that they use to extract strong- 
coupling interaction parameters. For a single layer stripe 
state, back- scattering interactions can be irrelevant only 
if the true ground state at v = 1/2 breaks particle-hole 
symmetry. This raises an interesting question. Could 
there be another class of as yet undetected phase tran- 
sitions that occur in the quantum Hall regime cither in 
the high N stripe state regime or for lower N? Broken 
particle-hole symmetry at v = 2 would imply a finite- 
temperature phase transition in the 2D Ising universal- 
ity class, for which the deviation of the Hall conductivity 
at v = 1/2 from e 2 /2h could be taken as the order pa- 
rameter. There is certainly no evidence for such a phase 
transition in experiment, although it might be washed 
out by disorder 36 even if it occurred. In any event, bro- 
ken particle-hole symmetry in the ground state would 
require that a phase transition occur between the high- 
temperature stripe state of Hartree-Fock theory that 
does not have broken particle-hole symmetry and a low- 
temperature stripe state in which particle-hole symmetry 
is broken and back-scattering is irrelevant. In light of the 
evidence that fluctuation corrections to the Hartree-Fock 
ground state are weak, we still believe that the simpler 
conclusion of our earlier work is more likely to be cor- 
rect, namely that particle- hole symmetry is not broken 
and that the smectic state is not stable. 

As we have emphasized several times the approach we 
have taken does not lead to standard coupled Luttingcr 
liquid behavior. In particular, the decays of fermion cor- 
relation functions at large distances, and of the tunneling 



density of states at small energies, are faster than power 
laws. This property occurs in our analysis because of 
broken translational symmetry in the stripe state which 
makes its energy functional invariant under a simulta- 
neous translation of all stripes. Barci et al. 24 , have ar- 
gued that this unusual property might signal a failure of 
the perturbative renormalization group transformation 
we have used, which rescales spatial coordinates along 
the stripe edges but not across them. Although we agree 
that our conclusions concerning the nature of the true 
ground state could in principle be altered if it were pos- 
sible to extend the perturbative RG analysis to higher 
order. Indeed this must happen when our analysis is ap- 
plied to low index Landau levels in which stripe states 
do not occur. We do not believe, however, that the un- 
usual correlation functions signal a greater likelihood of 
this eventuality than normally applies to lowest order 
perturbative RG calculations. In practice, the micro- 
scopic back-scattering amplitudes treated perturbatively 
are sufficiently weak in high Landau levels that our lowest 
order calculations seem likely to describe what happens 
down to the lowest temperatures available experimen- 
tally. In our view the approach we have taken should 
be trusted when experimental evidence suggests that the 
physics at low-energies is described by the stripe states 
of Hartree-Fock theory. 

We find here that the role of back-scattering interac- 
tions is quite different in the bilayer case compared to 
the single layer case. At very large layer separations, 
the single-layer case in which stripe state physics occurs 
down to very low temperatures for v ~ 1/2 is recovered. 
However, already for layer separations ~ 10£, we find 
that interlayer back-scattering interactions which drive 
the system toward a state with spontaneous interlayer 
phase coherence along the edges become important and 
lead to a state with a substantial charge gap. Our predic- 
tion of odd integer quantum Hall effects with anisotropic 
finite-temperature transport coefficients in surprisingly 
widely separated bilayer systems is an important result 
of this paper. This conclusion about the properties of 
spontaneously coherent stripe states in the absence of in- 
terlayer tunneling differs from that reached by Fertig and 
collaborators who, incorrectly in our view, ignore inter- 
edge coupling in considering the properties of coherent 
stripes. Interestingly, intra-layer back-scattering interac- 
tions that drive the system toward a Wigner crystal state 
are irrelevant in this regime. We conclude that stripe 
states are indeed stable in bilayer quantum Hall systems, 
unlike the single layer case, but not smectic metal states. 
It seems likely that for very small layer separations, back- 
scattering interactions will drive the system toward a uni- 
form charge density state with interlayer coherence, al- 
though our perturbative approach is not able to offer any 
substantial guidance in deciding this question. 

The study of stripe state physics in single-layer quan- 
tum Hall systems requires samples of exceptional quality, 
beyond that required for studies of fractional quantum 
Hall physics with lower index partially filled Landau lev- 
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els which can be studied at higher magnetic fields. It is 
still not possible to create bilayer quantum Hall systems 
with disorder that is as weak as that in single layer quan- 
tum Hall systems. Nevertheless recent samples appear to 
be of a quality that opens the physics of stripe states in 
bilayer systems up to experimental study. We expect on 
the basis of this work, and of previous theoretical work, 
that the physics will be rich, with much potential for 
surprises beyond the properties anticipated here. 
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